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We analyse the most recent experimental data regarding D° - D° mixing, allowing 
for CP violation. We focus on the dispersive part of the mixing amplitude, Mj^, 
which is sensitive to new physics contributions. We obtain a constraint on the 
mixing amplitude: Mf^l < 6.2 x 10~ n MeV at 95% C.L. . This constraint is weaker 
by a factor of about three than the one which is obtained when no CP violation is 
assumed. 



I. INTRODUCTION 

The ongoing searches for D° - D° mixing IB, ^ ||, ^, ^, S have not yet detected a 
signal of such mixing. Thus, the experimental data place an upper bound on the mixing 
amplitude. The value of this upper bound, however, depends on the assumptions one makes 
when analysing the experimental results. Specifically, the question of CP violation in D° 
- D° mixing has an important impact on the final answer. Most often, D° - D° mixing 
experiments are analysed assuming no CP violation. While this assumption is valid for the 
standard model, it does not hold for many new physics models. 1 (See, for example, the 
super symmetric models in [JT0|, [TT|.) Obviously, if the constraint on D° - D° mixing is to 
be used to test such new physics models, the experimental data should be interpreted in an 
appropriate framework ||12|| . We therefore present here the analysis of experimental results 
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1 In fact, since it has been recently suggested that the D° - D° mixing amplitude may be large even in the 
standard model fl , CP violation may be the most valuable clue for new physics in this system. 
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allowing for CP violation in mixing. 2 

The organization of this work is as follows: In section [TI| we present our formalism. We 
review the most recent experimental data in section [ITI| , and perform the analysis in section 
prV] . We conclude in section [V|. 



II. NOTATION AND FORMALISM 



We follow mostly the formalism of ref. [14]. The mass eigenstates are given by 



\D 1 p)=p\D°)±q 



(1) 



The mass and the width differences are parameterized as follows: 
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We define the D° and D° decay amplitudes by 

A f = (f\H d \D°) , Af = (f\ H a 



D° 



(4) 



and the complex observable A/ as 



P A f 



(5) 



In almost all models, CP violation in decay of the relevant modes can be safely neglected 
T3| , |Hf (for an exception, see |I5[1), leading to 



A f = A f - . 



(6) 



Now we can parametrize the effects of indirect CP violation in the relevant decay 
processes: The doubly- Cabibbo-suppressed (DCS) D° — > K + n~, the singly- Cabibbo- 
suppressed (SCS) D° — > K+K~, the Cabibbo-favoured (CF) D° — > R-7T+, and the 



2 We do not consider, however, CP violation in D° decays, which is absent in most new physics extensions 
0. 



three conjugate processes. We denote 3 
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where the and 5 are the weak phase and the strong phase, respectively, and 
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Next we define 



x' = x cos 5 + y sin 5 , 
y' = y cos 5 — x sin 6 . 

The rates of the DCS, SCS and CF decays are expanded for short times t ^ 1/T 



T[D (t)^K + n-] = e - rt \A K - 7T+ \ 2 

x R D + y/R^(l + A m )(y' cos (j) - x' sin <f>)Tt + 



1 + 2Am (y 2 + x 2 )(Tt) 2 
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T[D°(t) — > K-n + ] = e- rt \A 



K-TT-* 



X 



R D + VRd{1~ A m ){y' cos <p + x' sin <f))rt + -^(y 2 + x 2 )(Ttf 



T[D°{t) — > K + K~] = e - rt \A K+K -\ 2 x [1- (1 + A m )(y / cos0-x / sin0)rt] 



r [-D°(t) — ► R-n + ] = r D°(t) — ► # 



7T 



Several experiments measure the parameter yep, defined by 

)0 , 



Vcp 



t(D° — ► K'w 
t(D° — > K+K 1 



Note that A m in our definition is twice smaller than the A m used by CLEO Q . 
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with t being the measured lifetime fitted to a pure exponential decay rate for the specific 
modes [14| . If CP is a good symmetry in the relevant processes, this definition of yep 
corresponds to 

_ T(CP even) - T(CP odd) 
VCP = r(CP even) + T(CP odd) ' ( } 

since then the K + K~ state is an even CP state and the K~tc + state is an equal mixture of 
CP even and CP odd states. By fitting the decay rates in ( 15|) and fll3|) to exponents, and 
expanding for small A m we get fH}] : 

Vcp = V cos <p - A m x sin <p . (19) 

We are interested in the dispersive part of the mixing amplitude, M® 2 : Short distance 
contribution from new physics can affect M® 2 in a significant way. In terms of measurable 
quantities, | M® 2 1 is given by [16| 

2 ^4(A ! ^ ± ^W (20) 
1 121 16(1 -A m 2 " 



or, using eq. (||), 

\M D \ 2 - r 2 ^ + A ™ y2 (21) 
|Ml2 ' " r 4(l-^)- (21) 

III. EXPERIMENTAL DATA ON D° - W MIXING 

The neutral D system is studied by various experiments. First, the CLEO experiment @ 
measures the rates (0), (0): 

R D = (0.48 ± 0.13)% , 
y'cos^ = (-2.5^)%, 

x = (0.0 ±1.5)% , (22) 

9 4 - D OQ+0-63 
^ -^m — U-Z " J — 0.80 ' 

sin^ = 0.00 ±0.60 . 



The FOCUS experiment 0] provides a measurement of the ratio between the branching 
ratio of the DCS and CF decays. This measurement is consistent with CLEO data at the 
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TABLE I: Measurements of yep- 



Experiment Value 

FOCUS @] (3.42 ± 1.39 ± 0.74)% 
E791 § (0.8 ± 2.9 ± 1.0)% 

CLEO§ (-1.2 ±2.5 ±1.4)% 
BELLE § (-0.5 ± l-0±g;|)% 



BABAR § (1.4±1.0; 



-0.6> 
-0.7) 



70 



level of ~ 0.8a. However, as no direct measurement of the parameters is done, no stronger 
bounds on the parameters result. 

The value of yep is measured by the various experiments. Table | presents the various 
results. The world weighted average of yep is hence: 

y CP = (1.0 ± 0.7)% . (23) 



IV. INTERPRETATION OF THE EXPERIMENTAL DATA 

Our aim is to constrain the D° - D° mixing amplitude M^. First we combine (|l~2|) and 



(H) to get 

yep + A m sin (x' cos 5 + y' sin 5) = y' cos cos 5 — x' cos sin 5 . (24) 

The measured values of (|22"D and (^) can be used to constrain cos 5. Assuming first 4 A m = 
and also |sin0| w we find 

(1.0 ± 0.7)% = {-2.5t\i)% c os 5 - (0.0 ± 1.5)% sin 5 , (25) 

which implies a certain distribution for cos 5. Due to the sign difference between yep and 
y' and due to the relative smallness of x' it is expected that this distribution of cos S will 
be biased to negative values. By a full analysis, considering the measured values of A m and 
sin we can characterize the bias by stating the total confidence level value: 

cos5<0.7 (95%C.L.), (26) 



A similar procedure was followed in ref. EJJ|. 
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TABLE II: Comparison between mass and width difference parameters at 95% C.L. with different 



assumption on mixing parameters. 


Assuming cos 5 = l,cos <fi = 1 


No assumption 


\x\ < 2.9% 


\x\ < 6.3% 


-5.8% < y < 1.0% 


\y\ < 4.6% 



(and cos 5 < 0.0 at 68% C.L.). 

Since we have now a distribution for x', y' and cos 5, we may invert (|12|) to solve for x 
and y: 

x = x' cos S — y' sin 5 , (27) 
y = y' cos 5 + x' sin 5 . 



We note that the signs of x and y in ( |2"7| ) are not measured by current experimental results. 
Since the measured value for x' is distributed around zero the sign for y is determined by 
the sign of y' which, in turn, depends on the sign of cos 0. This sign is not provided by any 
measurement (all we know is that |cos0| ~ 1). Similarly, the sign of x is determined by the 
sign of both y' and sin 5, which are not measured. 

The resulting distributions for x and y are therefore in the form of two superimposed 
distributions for the two possible sign choices (denoted by the ± sign). We obtain: 

x ss (±2.8 ±2.5)% , (28) 
y P3 (±0.9 ±3.6)% . 



We note that these values are different from those quoted in [[H]] where it is assumed 
that 5 = <p = 0. When we consider the distribution of cos 5, the bound on x (and hence the 
bound on Amo) is weakened by a factor of about 2.2. The bound on y, however, (and hence 
the bound on AT) is strengthened. For comparison, table |TJ shows the 95% C.L. ranges for 
x and y in the two cases: One which assumes cos 5 = 1 and cos0 = 1, and one which takes 
the values mentioned. 

We evaluate now the D° - D° mixing amplitude. Taking the average decay width |T7| 



T D = (1.595 ±0.011) x 10~ 9 MeV , (29) 
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and using (|2lD , we obtain a distribution for M 12 which is maximal near zero: 

\M? 2 \ < 6.2 x 10" 11 MeV (95% C.L.) , (30) 

(and \Mg\ < 3.3 x KT 11 MeV at 68% C.L.). 

It is interesting to compare this value to the ones obtained by using some simplifying 
assumptions. First, assuming no CP violation in mixing, we set A m = but allow for 
S, ^ 0. We get 

\Mg\< 5.4 x Kr n MeV (95% C.L.) . (31) 
Second, we set A m = = and allow 5^0. We get 

\Mg \ < 4.0 x KT n MeV (95% C.L.) . (32) 
Third, we set 5 = 0, but allow A m , ^ 0. We get 

\Mg \ < 3.9 x 10" n MeV (95% C.L.) . (33) 
Last, we set A m = = 5 = and get 5 

\M&\< 2.3 x 10" n MeV (95% C.L.) . (34) 

This is the value which appears in fT7|| . Thus, allowing CP violation, the resulting constraint 
is about 2.7 weaker (i.e. larger) then the one which is obtained with the maximal set of 
assumptions. 

V. CONCLUSIONS. 

We interpret the most recent data from the experimental searches for D° - D° mixing. 
Allowing CP violation in mixing, we obtain the upper bound 

\Mg\ < 6.2 x 10- 11 MeV (95% C.L.) , (35) 

which is 2.7 times weaker than the naive calculation. 



Actually, it is enough to assume A m — 5 — since, in this case, the value of <f> affects only y, which docs 
not contribute to Mj§. 
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The actual upper bound for D° - D° mixing amplitude depends, therefore, on the model 
in question. Assuming that CP is conserved in D° - D° mixing, as is the case in the standard 
model, the bound is the one in (|32|) . (If, in addition, one is willing to assume that SU(3)- 
flavour symmetry holds in D decays, the bound is given by 034]).) For a more general model, 
with new sources of CP violation, eq. ( pop gives the present bound. Taking into account 
this weaker bound leads to modifications [jll|] compared to analyses that consider only the 



CP conserving bound |19|, EO, |2T 
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